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SUMMARY 


A 3ys tem for calculating the physical properties of 

soniUrStUr^wT 1 F ow wlth a * lal symmetry and super- 
rotational flow in a two-dimensional field was 

mined by use of the characteristics method. The 

fL-yifc WaS app ii ed bo the study of external and internal 
low for supersonic inlets with axial symmetry. For a 
circular conicai inlet the shock that occurred at the 
„ ° f, ' f became stronger as it approached the 

mS f • h ® inlet . and became a normal shock at the axis. 
*.1- . e § i0n for which strong shock occurred increased with 
increase oi the angle of internal cone at the lip of 

d f°J!. an lnlet with a central body the method 

. ^ hdr actenstics was applied to the design of an 

3hape that > theoretically, results in 
very efficient recorapression in the inlet: it was 
s own that if an effuser is connected with the diffuser 

he detemined? 1Utl0n Very Sma11 shock ' wave drag can 


INTRODUCTION 


The characteristics method for the determination of 

for e t!n ni - phe ? ome ? a waa first used by Prandtl and Busemann 
tor two-aimensional flow (references 1 and 2 ). For flow 
with axial symmetry Frankl (reference 3) used the method 
characteristics for determining the shape of a super- 
onic circular effuser with uniform exit velocity, and 
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B’errari (references Ij. end 5) independently used the char- 
acteristics method for determining supersonic phenomena 
for every type of boundary condition. Subsequently 
Guderley (reference 6 ) and Sauer (reference () transf ormed 
the system prooossd by Frankl and B'errart and obtained a 
different analytical solution of the problem. In all 
applications the hypothesis of potential flow was made; 
therefore the equation of potential flow was used. 

When shock waves that are not nlane (two-dimensional 
flow) or conical (flow with axial symmetry) occur in 
uniform flow the variation of entropy across the shock 
is not constant, and the flow behind the shoe x is no 
longer isentrooic and becomes rotational. Ij. the varia- 
tion of entropy is small, the effect of rotation of the 
flow is not important for determining the pressure distri- 
bution along a body and the theory of potential flow 
gives correct results. If the shock wave is strong 
has large curvature, however, the eifect of tne rotation 
becomes important and the flow must oe considered fio 
rotational. 


The method of characteristics can be extended to 
apply to rotational flow if, in place of the potential 
function for the differential equation of motion, the 
stream function considered by Grocco (reference c) is 
used. With the characteristics method for rotational 
flow a more exact determination can be made of the shape 
of the shock wave and the distributi on- of velocity and- 
pressure for phenomena in which the effect of rotation 
is important, as in the internal flow through supersonic 
inlets. The procedure of numerical calculation is similar 
and net much mere complicated Than that used for the c&jb 
of potential flow with axial symmetry. 


SYMBOLS 


p pressure 

p density 

3 entropy, mechanical units 

V veloci t y 

M Mach number 

V limiting velocity corresponding to adiabatic 

1 expansion to zero pressure 
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W 


c = 


u 


n 


X 

y 

P 

cp 

e 


l 


6 

0 

r 

R 

n 


speed of sound 

ratio of velocity to limiting velocity 


GO 


x-component of relative velocity 
y-component of relative velocity 

Cartesian coordinates 


{ 

K 


arc sm — 
M 


Mach angle 

\ » / 

angle between velocity V and x-axis, radians 


angle between tangent to shock and direction of 
velocity of flow in front of shock 

deviation of direction of velocity across shock 
wave 

angle of polar coordinate in conical field 
potential function 

stream function for rotational flow (see equa- 
tion ( 11 ) ) 

ratio of specific heat at constant pressure and 
constant volume 

ga3 constant 

normal to streamline 


j - 3 ^- n _P tan P sin cp 

cos (cp + p) 

_ sin p tan p sin cp 

cos (cp - p) 

II , L, K, and N defined by equation (lip) 
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£l = X A " X B 
S 0 = X C " X A 



Wg 


g = 


1 + 



v* 


sin p A 


COS 

(P A 

+ 1>A> 


sin 

Pb 

cos 

(Vr 

" Pb) 

tan 

Pa ' 

+• h tan 

s 3 ' 

- s A 

1 


YR e + gf 


Subscripts : 

0 chamber condition (zero-velocity adiabatic trans- 

formation) 

A points of first family 

B points of second family 

C quantities in the points calculated from A 

and B 

x derivative with respect to x 

y derivative with respect to y 

a ahead of shock 

b behind s ho ck 
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cpp value corresponding to valuo of cp at point C n 

n 

C n at point Cn 

CHARACTERISTICS METHOD FOR SUPERSONIC POTENTIAL 
FLOW WITH AXIAL SYMMETRY 


The differential equation for potential flow with 
axial symmetry (reference )+) is 

= o (1 

a 2 / ox 2 \ a 2 /dy 2 a 2 dxdy J 

In supersonic flow some lines can he individuated 
(characteristic lines) that divide the flow into two 
regions for which the values of 0, 0 X , and 0j along 

the line are different. For every point of the flow two 
characteristic lines can be determined; every line is 
inclined at the Mach angle with respect to the direction 
of the velocity at the point, and therefore the character- 
istic lines can be divided into two families on the basis 
of the sign of the angle of the characteristic line with 
respect to the direction of the velocity. A family that 
is usually called the first family Is defined by the equa- 
tion 



dy 

dx 


tan ( p + m ) 


( 2 ) 


and the other family (second family) i3 defined by 


dy 

dx 


tan ( cp - (3 ) 


(3 


The variation of the quantities that define the 
velocity ( and V) along a characteristic line is 
given by the following equations from reference 5: 
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For the first family, 

~ - dcp tan p - l-y- =0 (1+) 

and for the second family. 


Y + dcp tan p - m— =0 (5) 

where l and m are trigonometric expressions defined 
as 


sin p sin cp tan p 
cos ( p + cp ) 


( 6 ) 


_ sin p sin co tan p 
cos (cp - p) 


If the direction of the velocity and the Mach number 
at two points near each other (points A and B in fig. 1) 
are known, the direction of the velocity and the Mach 
number at a point C given by the intersection of the two 
characteristic lines of different family starting from 
points A and B can be determined. Because the distances BC 
and AC are small, all the coefficients of equations ( q) 
and (5) can be considered constant and coincident to the 
corresponding values at points B and A, and. the tangents 
to the characteristic lines at points A and B can be 
substituted for the characteristic lines from point A to 
point C and from point B to point C. In this case the 
lines AC and BC are straight lines. The line BC is 
inclined at an angle cp - p with respect to the x-axis 
and the line AC is inclined at an angle cp + p with 
respect to the x-axis. Equation (ij.) can be applied for 
the line AC where 


I 
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dX = Xn “ 

y = 

V = V A 

P = Pa 

i = i A 


and equation ( 5 ) can be applied for the line BC where 

dx = Xq - x B 


y = YB 
V = V B 

P = Pb 


m = m B 


For practical use equations (k) and ( 5 ) are combined and 
transformed by means of the ratio W of the velocity V 
to the limiting velocity Vj. This ratio is defined as 



= 1 + — - — sin 2 13 (7) 

Y - 1 


and the following equations are obtained: 
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/ \ W A 

dcp A f tan p B + r— tan M = 1 “ ~ - (<Pa “ ^ 5 ) tan 


np 


YB ( 


x c - x b) - ( x c - x a) ^7 v 


Z.A WjA 
A v; B 


(8) 


dW 


A 


y A 


= tan p A dcp A + 


(Xg - X A )l A 


( 9 ) 


Vvg — W A + dW A 


<?C = <Pa + dQ5 A 


10) 


With the method of characteristics (reference J+) it 
can be shown that if a deviation of a streamline which 
wets the body occurs, the phenomena on the corner are 
regulated by* the same laws that regulate the two-dimensional 
flow; therefore, the tangents of every characteristic line 
starting from the corner are known. If the initial flow 
conditions are known, the step-by-step calculation of all 
the physical properties of the flow in the entire field 
is permitted, particular ly the calculation of the shape 
of the shock wave and the pressure distribution along any 
body of revolution with axial symmetrical flow In cases 
in which the hypothesis of potential flow Is correct. 


CHARACTERISTICS METHOD FOR SUPERSONIC ROTATIONAL 
FLOW WITH AXIAL SYMMETRY 


Supersonic perfect flow is rotational when the flow 
is preceded by a shock wave and when the variation of 
entropy across the shock wave changes from point to point 
behind' the shock. In this case the transformation of the 
fluid along every streamline is isentropic until another 
shock wave occurs in the fluid (reference 8). 


I 
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If a stream function \j/ is assumed to be defined 
by the following equations from reference 8: 


1 , 

\{/ y = yu( 1 - w2) Y_1 


1 

= -yv( l - w2) Y_1 

the equation 




> 


(id 


f U /) 



Y-l 

y(i - w 2 ) 


( 12 ) 


is a function of only the stream function (reference 8) 

and therefore from the equations of state, continuity, 
energy, and steady motion the following equation can be 
obtained: 




- y 2 ( l - w2) 



fty) =0 



Equation (13) is a Monge-Amoere equation, and if 

,2 


H = 1 - 


u £ 


.2 


L = 1 - - 

C^ 


K 


N 


uv 


= i. y 2 a . w 2 ) r-i(^.^ 


f(^) (111) 
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two characteristic families with the following equations 
can be obtained: 


Por the first family, 



and for the second family. 


£y _ K /k£ _ L 
dx: H \ H 



(15) 

(16) 


(17) 


(18) 


If n is the normal to the streamline, equation (11) 
yields 

o 

•^^2 + - y2\fy 2 ( 1 _ Vi/2 ) ^ 


and 


— = grad2 \|/ ( 19 ) 

w 


curl V x V 


a 


2 


— grad s 

yr 


1 ds 
YR dn 


(20) 


because 


s is constant along every streamline; therefore 
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f(*) = 


Y - 1 1 ds 
2 Y R dn 


yW( 1 - W 2 ) Y " 1 


/ 


21 ) 


and the following expressions can be obtained: 




v xx 


yv 


x 


(1 - . «2(i - wpl- 


ax 

-Y* 2 


v 


+ SJ2L(1 - *2) ^ *, 

Y - 1 3 


4- = yuxd - + Si ud - vdpmr 

y x dx 

- i+2 

- H2L,! - *2) ^ w 

r - i x 


( 22 ) 


,,T hon equations (21) and (22) are substituted in equations 
(15) and (±7) and (l6) and (l8) and the Mach angle and 
the velocity are expressed in polar coordinates, equations 
( 13 ) and ( 17 ) become 

^ = tan ( P + cp ) (first family) ( 25 ) 

dx 

= tan (cp - p) (second family) (2li) 

dX 

and equations (l6) and (l8) become 
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dW L a . dx 

t an p dqp - 

W 


dx as sin* ,3 


y 7 dn yn oosTp +" c ) 


= 0 (first family) (2p) 


dW 


+ tan (;• do - — m - — ~ ^ — *-r- = 0 (second family) (26) 

y dn yu cos ( cp ~ P/ 


donations (23) and (2l+) are identical to equations 
for potential flow (2) and (3), and equations (2S) and 
(26) rare similar to equations for potential flow \b.) and (S' > 
differing onlv by the terms that contain ds. Equations 
(25) to (?. 6 ) permit a step-by-step calculation of the 
entropy, intensity of velocity, and direction of the 
velocity if the initial and boundary conditions are known. 

If all the physical properties are known for two points 
A and B and if the two points are close to each other, 
the tangents to the characteristic lines at the points A 
and B can be substituted for the characteristic iines^ 
wi th close apo ro;;l mat! on. In tnis way a point C can os 
determined as the intersection of the second characteristic 
line of point B with the first characteristic line of 
point A (fig. 1) , because cp and p are known for the 
points A and B. 

For the characteristic lines of the first family, 
equation (23) gives the variation of cp and W from 
point A to point C, end all the coefficients are known 
and correspond to the coefficients for point A; only the 

term —2. is unknown. From equation ( 26 ) the variation 

dn 

of cp and W from point B to point C can be determined, 

as 


dn 


are known and equal 
d: 


and all other terms, except 
tc the values for point B. The term 
from the value of the entropy for points A and 3 . 


can be determined 


dn 


1 
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From figure 2 , in equation (25) the term 


ds 

dn 


csn oe 


written as 


ds ( a C " 3 a ) coa (Pa + <Pa) 


dn 


(x c - x A ) sin p A 


( 27 ) 


and in equation ( 2 o) can be written as 


d 3 ( 3 C - S3) cos (cpb - Pb) 


dn 


(x c - xg) sin p B 



If the points A and B are close to each other and the 
variation of entropy is gradual, equation (25) can he 
written as 


ds 

dn 


3 B 


(*C “ X A) 


sin p A 


cos (Pa + <P A ) 


3 A 


+ (-C - *b) 


s m 


P'B 


(29) 


cos (cp 3 - p B ) 


and - — can be considered eaual in equations ( 23 ) and (2 
dn 

d s 

Tn this case — is known; therefore cp and W can be 
dn 

calculated for ooint €. For practical calculations equa- 
tions (25) and (26) should be transformed Into two equa- 
tions (equations (Jl) and (32)) each of which contains 
only one of the unknown terms dw and dcp. 


VO 
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For simplicity, let 

£>1 ~ x a “ X E 
So - X C " X A 
Si 


= 1 + 


"O 


A CD 


<?A ’ s°B 


W, 


h = 


e = 





s in 

Pa 


COS 

(Pa 

+ 

<pa> 


sin 

Pb 


CO s 

(<PB 

- 

Pb) 

= tan 

Pb ' 

+ h tan 

S B 

- S A 


i 






( 50 ) 


Y R e + S-f 

Then the following equations can )?e obtained- 


0 „ , Sc^bS So hJ A 

r drn A = 1 - h - Aep tan p B +A(gf sin 2 f B -he sin^p B ) + “ (31) 


Yb y A 


d'vVA. 


= tan (3^ dcpA + Ao ain^p^ + 


W 


A 


SO* A 
Ya 


( 32 ) 


i 
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3 C ~ S A + AeyR 


<?C = <?a + dc P; 


W q — + dw^ 


sin^ 


- r - 1 



(35) 

(3k) 

(35) 

(36) 


^0 m B® So hl A 

and become 


In equation (31) the terms 

y B y A 

very imnortant as y approaches zero, near the axis of 
the body; therefore, in this region the distance between 
the points considered must be reduced to obtain sufficient 
accuracy . 


CHARACTERISTICS METHOD FOR TWO-DIMENSIONAL 
SUPERSONIC ROTATIONAL FLOW' 


In two-dimensional flow an equation similar to equa- 
tion ( 13 ) can be obtained if a special stream function 
defined by the following equation is assumed (reference 8); 


i'y = u(i - w*) 


- 1V2>Y-1 


'!/ x = -v(l - ^ 2 ) Y-1 


(37) 


In this case the equation of motion (equation (13)) becomes 
( reference 8 ) 



( 3 2 ) 
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Equation (38), like equation ( 13 ), is a Monge-Ampere 
equation and permits the determination of two equations 
that define the characteristic lines and two equations 
that give the variation of the velocity along the char- 
acteristic lines. The equation with transformations 
analogous to the case of three-dimensional flow can be 
written in the following form: 


dy 

— = tan ( p + cp) 

ClX 

dy 

— — = tan (o - p) 
dx 


(first family) 
(second family) 


( 39 ) 

(ko) 


diV , , dx ds sin3p 

tan p dcp — 1 = 

IV yR dn cos (p + cp) 


0 (first family) (I|l) 


d'.V dx ds s in* p , 

— + tan 3 dcp- : — = ° (second iamily) (R2) 

W vR dn cos (cp - p) 

Equations for cp and aW similar to equations ( 31 ) 
and (32) can he obtained from equations (39) to (^-2) 

using equation ( 29 ) for the expression — as follows: 

dn 


r dcp A = 1 - h - Acp tan pT 


dV. A 


A 


sm‘ 


■Pc 


3 + n(gf sin 2 p B - he 

sin 2 p A ) (Il3 ) 

p A dcp + Ae sin 2 p A 

(U+) 

= s A + AeqR 

a, -5) 

= cp A + dcp A 

( 16 ) 

= W A + d, A 

( 1 - 7 ) 

: ( A- . A 

(L8) 

2 Vr 2 ) 


l 
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DETERMINATION OF SHAPE OF SHOOK AND PRESSURE 
DISTRIBUTION ALONG A BODY 


The physical properties of supersonic flow pasu a 
hod 7/ of revolution in axially symmetrical flow can e 
determined step by stop by tno use of equations ^ 

to (j6). The system of calculation is similar to tnuo 
for potential flow. If the Body begins with. a poin. 
(reference 5 ), a cone tangent to the body (fig* 2 ) can 
he substituted for the front part of the body. If point a 
is the point at which the body can no longer be considered 
coincident with the cone, the velocity at point A is knov/n 
from the cone calculations (references 9 and 10 ); tnero- 
fore, the shape of the characteristic line AB of the xirst 
family can be designed, because <p and p corresponding 
to different angles 6 are known from the cone calcula- 
tions. At point A the body turns through an angle Acp 
and the flow undergoes a transformation that is deter- 
mi nab la bv the laws of two-dimensional flow, ihe velocity 
and direction of the flow after the deviation Acp there- 
fore, and the tangents to the new characteristic lines 01 
the first family at point A can be designed. At a poan- , 
near point A, the intensity and direction of^the^ velocity 


m own : cons e que nt 1 y , t h< 


are 

line of the second 
the point Cq can he 

properties can be c 


can 

d , 


tial flow (equation 


fro: 


it of the 


body 


amily 
loc 

Iculated 
( 8 ) and 
is conical 


tangent to tlio characteristic 


point v^ n 
In order 
wave , th 


Cv, is determined from 


be designed at point B and 
At point C 1 the physical 

h the equations of poten- 
be cause the shock in 
. with the same system the 
the rsoint B v aV,r ‘ rtV 


vv 1 

/ r 

\ y 


)) 


to determine the flo. 

) equation across the 
1 


point 
the 


r 01 
shock 


the 
If 

the 

6 the deviation of 
shock ’wave , and the 
conditions ahead of 
the shock equations 
(reference 11): 


, n on the shock wave . 
point C n+1 on the shock 
and the equation for 

characteristic lines of tho first family must oe used. 

€ is the angle between the tangent to the shock and 
direction of velocity of flow in front of the shoex, 

the direction of velocity across the 
subscripts a and b denote the 
and behind the shock, respectively, 
can be written in the following form 


tan e 


2 


tan (e - 6) y + 1 


Mb 2 


sm^ 


2f e 


- 6 ) 


- 1 


(^- 9 ) 
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MIA Id Ao . 


Y + 1 


M q 2 


tan 5 


1 3^2 sin^e - 1 J 


- 1 ) tan e 


(50) 


5 b 2r / n n Y “ 1 

— = r ( M a ^ sinM - TT~ 

o- v + 1 V ° AY 


-a 


Y + 1 \ 


(5D 


a 2 


P )3 Y + 1 V I.! a 2 sin 2 « 


1 + x^l 


(52) 


s b “ s a 


r / p b . p a\ 

—“7 lo Se7T + Y los e7“ ; 
Y - 1 V p a 


,p b J 


(53) 


— = i + — - — — 

■,2 Y — 1 5 f 2 

“ 1 i'JV', 


(54) 


Equations ( 5 1 ) to (55) show that the general condi- 
tion for obtaining potential flow behind the shock when 
the flow in front of the chock is potential is that M a sin e 

must be constant. Since the value of I’ a in front of 
the shock is known from equations (1+9) to (54)> the 
d( 3b - s a ) dv .r 

values of and — can be determined as 

d5 d5 

functions of b and therefore the equation for the char- 
acteristic lines of the first family (equation ( 2p ) ) 
becomes 


VVcp C n " ' l ’ C n ^ dco 


OX 


sm 


3p c 


+ 


W, 


dfi 1 tin 


- tan (3c n dc o Iq 


n 


'n 


n 




YR 


fCn - S C 


n 


d ( s-^ - s a ) 


d.6 


dco 


= 0 


(55) 



NACA TN No. 1135 


19 


where W n is the velocity behind the shock, 3 ffln is 
s c n a 0 \ vc, n 

(i W Cl ( S “ S ) 

the entropy, — - and are coefficients that 

d6 d6 


correspond to a deviation across the shock for which the 
velocity assumes the direction of cpQ n , ’^Cn ^he 

velocity, and SQ n is the entropy at the point C n . Vi/hen 


the value of d<p has been determined, the deviation 
across the shock 6 and the corresponding values of W 
and e can be determined at point C n+ p. 


In order to determine the velocity on the body at a 
point A 1 , equation (26) is used. At point Ap the value 
of cp is known because the direction of the flow is 
tangent to the body and, therefore, dcp is known; 
ds 

— is also known because in equation (27) the value 
dn 

of sq corresponds to the value of s at the point Ai 
and is equal to the value at point A, which is known from 
the cone calculation. Prom the value of the pressure 

on the body relative to the pressure p^ T (pressure for 

zero velocity from isentropic transformation from the 
conditions behind the shock) can be determined as 


Y 


z~T = (1 - w2 ) 

Po 


r-i 


(56) 


On the lip of the nose of an open-nose body a shock 
wave that is a two-dimensional shock occurs; therefore, 
the tangent to the shock on the lip is known and the 
pressure and velocity behind the shock are also known 
(equations (1;9) to ( 5 I 1 ) ) • In figure A the line A3 and 
the velocity behind the shock at point B are known. With 
equation ( 26 ) the velocity at point C can be calculated 
with the system just described, and from the point B the 
point D along the characteristic line of the first family 
can be calculated by equation ( 55 )* Prom point D the 
point E can be calculated and a point P can be inter- 
polated, which permits the determination of points G and H. 
Because the curvature of the shock near the lip is large, 
a point C that is very near point A must be selected so 
that correct interpolated values may be obtained. The 
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point E can be recalculated from point G, and the points I 
and M can be recalculated from point F in order to have 
a second approximation. 

For the case of two-dimensional flow the procedure 
is the same as that for three-dimensional flow; equa- 
tions (!p5) to (Ipj) are used instead of the corresponding 
equations (equations (pi) to (5b)). 

EXAMPLES OF APPLICATION OF CHARACTERISTICS METHOD 
FOR ROTATIONAL FLOW 

Determination of Shock Shape and Pressure Distribution 
along the External Surface of a Slender 
Open- Nose Body of Revolution 


Theoretical and experimental calculations were made 
to determine the shock shane and pressure distribution 
along the external surface of a slender open-nose body of 
revolution. The body considered is the same body for. 
which calculations of* the external pressure distribution 
were made by Brown and Parker of the Langley Memorial. 
Aeronautical Laboratory by use of the snail-disturbance 
theory. The calculations were made for a free— stream 
Mach" number of 1.525, for* which a schlieren photograph 
taken during tests was aval lab X 
calculated, and test results . 


lor comparison of the 


In order tc determine the importance of rotation of 
the flow, calculations for the front part of the body were 
also made with the potential-flow characteristic equations. 
The results of these calculations differed only Slightly 
from those obtained from the characteristics method for 
rotational flow because the curvature of the shock was 
very small. The calculations were begun with the deter- 
mination of the two-dimensional shock: on the lip of the 
nose. The practical system of the calculations i3 shewn 
in figure Ij., and in figure 5 the calculated streamlines 
and shock-wave shape are compared with the shock-wave 
shape obtained from test results. In figure 6 the pres- 
sure distribution calculated by the characteristics method 
is compared with the pressure dis trioution determined b y 


i 
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the sraall-cti sturbance theory. The small-disturbance 
theory undervalue s the Increase in pressure 
that occurs through the shock, but the differences in the 
results obtained by this method and those obtained by the 
characteristics method are small. 


Determination of Shock Shape, Streamlines, and Pressure 
Distribution along the Internal Surface of a 
Slender Open-Nose Body of Revolution 

Three slender open-nose bodies with different conical 
inlet angles are considered (figs. 7 to 12) and the super- 
sonic part of the internal flow is studied for a free- 
stream Mach number of 1.6. For this type of body the 
internal shock produced at the lip of the inlet has a very 
large curvature and the effect of rotation is therefore 
very important. The calculations are extended to the 
region in which the Mach number is 1.0. The hypothesis 
is made that subsonic boundary conditions and stability 
considerations permit a subsonic flow at the end of the 
supersonic flow such as results from the calculations. 

The results of the calculations show that at the axis of 
the inlet a normal shock always occurs and that the region 
in which a strong shock occurs (with subsonic velocity 
behind the shock) increases in extension with the increase 
in angle of the internal cone. When the cone angle 
approaches zero and the shock is a Mach wave, a complete 
reflection occurs at the axis of the inlet and the exten- 
sion of the strong shock is zero. 

For large internal cone angles (figs. 9 to 12) the 
shock is a simple shock that becomes normal in the central 
part of the body of revolution. After the shock the 
compression continues but the characteristic lines cannot 
form an envelope. For small internal cone angles com- 
pression occurs gradually and an envelope of Mach lines 
occurs. The shock therefore reflects from the central 
part and another shock is generated. The form of the 
reflected shock is shown in figure 8. The ratio of the 
diameter of the region in which a strong shock occurs to 
the diameter of the inlet as a function of the internal 
cone angle is shown for M = 1.6 in figure lp . The 
results are interesting for the practical design of 
supersonic inlets of slender shape because they show that 
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for large cone angles the central part of the body of 
revolution, in which the compression Is not very efficient, 
is large. 


Determination of Physical Properties of the Internal plow 
through an Inlet with a Central Body 

Por a Mach number of 1.6 an analysis ol the shape .of 
an inlet with a central body was made to aid in obtaining 
high efficiency. Theoretically a supersonic diffuser with 
or 'without a central body and having no shock losses or^ 
shock drag can be obtained (reference 12 ; ; out for pracuical 
use it is convenient to accept small shock drag in order 
to avoid large friction drag. The inlet considered 
(fig. lip) has a 10° central cone. The deviation across 
the conical shock is 1+8 ? . The shock is reflected by a 
cylinder that forms the external part of the inlet. The 
reflected shock produces rotational flow behind the shock 
and the deviation across this shocu on the cylinder is 1 
and on the central body, 2°Lp3 ' . If the central body, 
behind the shock has the same direction as the velocity, 
the shock will not be reflected and isenorooic compression 
can be obtained behind the shock (fig. l4)» design 

of the central body, therefore, was determined from the 
calculation of the’ corresponding stream 
tion of the velocity along 
value of the exit 
condition and the 
velocity at every o 
errors stream tubes 
basis of 
velocity 
check on 
The Mach 
fixed at 
f rom the 
instability. 


_ine. The varia- 
te external cylinder ana the 
velocity were fixed, and from this 
condition dependent on the shock, the 


int was calculated. In order to avoid 
were designed that permittee, on the 
the ratio of the area in tnc region ol unilorm 
to the area at the eno of the stream uube, a 
the precision of the numerical calculations, 
number in the minimum section of the inlet was 


value larger than 1.0 so that 


subsonic 


The 


cart of the flow would 
value chosen was I. 09 . 


disturbances 
not cause 


If an effuser is connected witn the aiffuser,.a body 
of revolution with very low shock drag can be obtained 
(fig. Ip). The only pressure losses are the losses . across 
the two shocks, which are very small; but for practical 
applications the friction losses are larger tnan for the 
internal body alone. A balance of the pressure losses 


I 
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and friction losses must therefore "be made in order to 
examine the possibility of practical use of this arrange 
rrtent . 


CONCLUSIONS 


A system for calculating the physical properties of 
supersonic rotational flow with axial symmetry and super- 
sonic rotational flow in a two-dimensional field was 
determined by use of the characteristics metuod. practical 
use of the system is based on a step-by-step procedure, 
which requires long numerical calculations; but the 
calculations for three-dimensional flow are of the same 
type as for potential flow and, therefore, can be used 
for the practical problems in which rotation is important . 
Some applications were made to determine the external and 
interna] flow on bodies of revolution with axial symmetry, 
and the following conclusions were indicated: 

1. The effect of rotation is not very important if 
the variation of entropy is small but is important in the 
study of internal flow, for which the variation of entropy 
is usually large. 

2. Wien the inlet .is a circular ■ conical channel, a 
shock is produced at the lip of the inlet that becomes 
stronger when the shock approaches the axis of the inlet 
and becomes a normal shock at the axis. The region for. 
which a strong shock occurs (with subsonic velocity behind 
the shock) increases with the increase of the angle of 
internal cone . 


3. If an inlet with, a central body is considered, 
the method of characteristics permits the cesign of an 
internal- channel shape that, theoretically,, results in 
very efficient re compression in the inlet; and, if an 
effuser is connected with the diffuser, a body of 


revolution with very small shock-wave 
mined . 


can bo detor- 
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APPLICATION OF THE METHOD OF CHARACTERISTICS 
TO SUPERSONIC ROTATIONAL FLOW 
By Antonio Ferri 
September 1946 


Page 12, equation ( 25 ): The last term on the left-hand side of 

equation ( 25 ) should be preceded by a plus sign instead of a 
minus sign. 

Page 13, 

line 1 : "Figure 2 " should be "figure 1 . " 

equation (28): (sg - eg) should be used instead of (sc - sp) 

in the numerator on the right-hand side of equation ( 28 ). 
line 2 under equation (28): "Equation ( 25 )" should be 

equation ( 27 )." 

Page 14, 

equation ( 31 ) • The fourth term on the right-hand side of 
equation ( 31 ) should read as follows: 

A(gf sin 2 3B + be sin^P^) 

equation ( 32 ): The second term on the right-hand side of 

equation ( 32 ) should be preceded by a minus sign instead of 
a plus sign. 

Page 16, 

equation (4l): The third term on the left-hand side of equation (4l) 

should be preceded by a plus sign instead of a minus sign, 
equation (43): The expression in parentheses in the last term of 

equation ( 43 ) should contain a plus sign instead of a minus sign; 
thus 


A(gf sin^PB + he sin 2 p^) 

equation (44): The last term on the right-hand side of equation (44) 

should be preceded by a minus sign instead of a plus sign. 

Page 18 , equation (55): The factor before the brackets in the last 

term on the left-hand side of equation ( 55 ) should be preceded by a 
plus sign instead of a minus sign and should contain sin 2 instead 
of sin^j thus 

sin 2 Pc n 
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